On the Bose— Einstein distribution 
and Bose condensation 
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Abstract. For a system of identical Bose particles sitting on in- 
teger energy levels, we give sharp estimates for the convergence of 
the sequence of occupation numbers to the Bose-Einstein distri- 
bution and for the Bose condensation effect. 
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1. Introduction 

The Bose-Einstein distribution has been studied by physicists (here 
we only mention the well-established textbooks by Landau and Lif- 
shits [1] and Kvasnikov [2]) as well as mathematicians (e.g., see Ver- 
shik's papers [3,4], where further references can be found). It is well 
known that the vector of occupation numbers tends, in a certain sense, 
to the Bose-Einstein distribution as the total energy of the system and 
the number of particles tend to infinity. It is, however, of interest to es- 
tablish sharp estimates for this convergence, including the case of Bose 
condensation. This has been done in our papers [5-9]. The present 
exposition is mainly based on these papers but contains a number of 
improvements and has the important advantage of being largely self- 
contained (at least as far as Bose condensation itself and distributions 
with variable number of particles are concerned). Moreover, in con- 
trast to other mathematical treatments of the subject, it does not use 
any but very elementary mathematical tools. (For example, we avoid 



2000 Mathematics Subject Classification. 60C05 (Primary) 82B30 (Secondary). 

1 



2 



V.P. MASLOV AND V.E. NAZAIKINSKII 



resorting to number-theoretic results like the Meinardus theorem [10] 
or methods of analytic number theory [11].) 

Before proceeding to the results themselves, let us recall what Bose 
condensation is by using an elementary combinatorial model, namely, 
that of balls distributed over boxes. 

Suppose that there is a sequence of boxes Uj, j = 0, 1, 2, ... , and 
each box Uj is divided into qj compartments. We take N identical balls 
and put them into the boxes at random observing the only condition 
that 

oo 

(1) J]jiV,<M, 

j=0 

where Nj is the number of balls in the box Uj and M is a positive 
integer specified in advance. As an outcome, we obtain a sequence of 
nonnegative integers Nj, j — 0,1,2, ... , such that 

oo 

(2) E^:'-^^ 

j=Q 

and condition (1) is satisfied. It is easily seen that, given M and N, 
there are finitely many such sequences. Suppose that all allocations of 
balls to compartments are equiprobable. Since the number of ways to 
distribute Nj indistinguishable balls over compartments is equal to 



Nj J Nj\{q, - 1)! ' 
it follows that each sequence {Nj} can be realized in 

(3) m{N,}) = fi('''^''-') 

j=0 ^ 3 J 

ways, and the probability of this sequence is equal ioWiXNjY) divided 
by the sum of expressions similar to (3) over all sequences of nonnega- 
tive integers satisfying the constraints (1) and (2). This makes the set 
of all such sequences a probability space; the corresponding probability 
measure will be denoted by ^m,n- The positive integers qj are called 
the multiplicities. We will assume that 

(4) q, = Qf-' + o{f-'), j-oo, 

where o? > 1 is a given parameter (which we refer to as dimension) and 
(5 > 1 is a positive constant. 

What happens as M, N — > oo? It turns out that the so-called 
condensation phenomenon can occur: if N tends to infinity too rapidly, 
namely, if it exceeds some threshold N = N[M), then a majority of 
the excessive N — N balls end up landing in the box Uq; more precisely, 
with probability asymptotically equal to 1, the number of balls in Uq 
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is close to N — N (and accordingly, the total number of balls in all the 
other boxes is close to A^, now matter how large itself is). Let us 
state the corresponding assertion (which is a special case of Corollary 4 
and the subsequent argument in Sec. 2.2.2). 

Theorem 1. Define N = N{M) by the formula 

oo 

where b is the unique positive root of the equation 

oo 

and set 

' (N In Ny/^x(N) ifd>2, 

N^^"^ In Nx(N) ifl<d<2, 

where x(.x), x > 0, is an arbitrary positive function arbitrarily slowly 
tending to infinity as x ^ oo. There exist constants Cm, independent 
of M and N, such that if N > N, then 

Pm,n{\No -{N-N)\>A)< CrrJ^'"", m = 1, 2, . . . . 

Nothing of this sort happens if < A^. In this case, the limit distri- 
bution as M, A^ — > oo is the Bosc-Einstcin distribution with parameters 
(3,11 > (see formulas (22) and (23) below), and no condensation on 
the zero level occurs. 

It is not hard to write out an asymptotic formula for A^. To this end, 
one substitutes (4) into (5) and (6) and applies the Euler-Maclaurin 
formula to the resulting series so as to transform them into integrals. 
The result (see formula (19) in Sec. 2.1.2) is that 

(7) N = C{d)M^Q^{l + o{l)), 

where C{d) is a constant^ depending only on the dimension d. 

Suppose that, for given M and N > N, we wish to avoid Bose 
condensation. How can we do that? 

One way would be to increase A" so as to ensure that N > N. 
To this end, let us partly "Boltzmannize" the system, i.e., make the 
balls partly distinguishable. More precisely, suppose that the model 
is basically the same, but we are additionally allowed to paint each of 
the N balls at random into one of K distinct colors. Now that we can 



^The explicit expression is 

(r(d + i)c(rf + i))^ 

where T{x) is the gamma function and ({x) is the Euler zeta function. 
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distinguish between balls of different colors (but balls of a same color 
are still indistinguishable), the Bose condensation threshold should 
change. 

Let us compute how exactly it changes. To make the computa- 
tion, instead of painting the balls, we mentally divide each of the qj 
compartments in the jth box into K sub-compartments and put the 
uncolored balls there (with the understanding that putting a ball into 
the kih. sub-compartment is equivalent to painting the ball into the kth. 
color). Now there are Kqj compartments in the jth box, and we see 
that, all in all, the introduction of K colors has the only effect that all 
multiplicities qj are multiplied by K. 

Let us apply Theorem 1 (with qj replaced by the new multiplicities 
qj — Kqj). Formula (19) gives an asymptotic expression for the new 
threshold, which we denote by N . All we have to do is to replace Q by 
KQ in formula (7); then we obtain 

N = C{d)M^{KQ)^{l + o(l)) = K^N{1 + o(l)). 
Thus, the introduction of K distinct colors has raised N by the 



In this section, we state our main results. All proofs are given in 
Sec. 3. The simplest physical model to imagine behind our mathemat- 
ical constructions is that of a system of identical Bose particles sitting 
on integer energy levels. 

2.1. System with a variable number of particles. Wc start 
our analysis by considering the "photonic" case, where the number of 
particles in the system is not fixed and only a constraint on the overall 
system energy is given. 

2.1.1. Definition of the system. Let M > be an integer. We 
denote the set of all sequences {Nj} = {A^j}^i of nonnegative integers 
satisfying the condition 



by flM- Note that all such sequences are finitely supported^ and Qm is 
finite. 

Next, let positive real numbers qj > 0, j = 1, 2, . . . , be given. We 
introduce the probability space Xm = {^m, ^m, Pm), where J-'m = 2^'^" 
is the powerset of VIm (as is customary with discrete probability spaces) 



factor Kd+i . 



2. Main Results 



oo 



(8) 




^That is, all but finitely many Nj are zero. Indeed, Nj = for j > M. 
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and the probability Pm is defined as follows. We assign the weight 

j=l \ J / 

where (^J is the generalized binomial coefficient,"^ to each element 
{N.j} G CIm and the weight 

(10) w{A)= ""({^.» 

{Nj}eA 

to each subset A C 0,m and set^ 

(11) p^(^)^J^, AGflM. 

The numbers qj will be referred to as multiplicities.^ We assume 
that they have the asymptotics^ 

(12) qj = Qf-\1 + o{l)) asj^oo 

for some real constants d > 1 and Q > 0. In particular, there exists 
constants Bi,B2 > such that^ 

(13) B,f-'<qj<B2f-\ j = l,2,.... 

2.1.2. Limit distribution. It turns out that if M is large, then a 
randomly chosen sequence {Nj} G Qm is with high probability close, 
in the sense described in Theorem 2 below, to the nonrandom sequence 
{Nj} defined as follows. For M > 0, the equation 



OD 



(14) 5:^^=^ 



"'Recall that (^^J = nj=c[ empty product (for n = 0) being by definition 

equal to 1. Since {Nj} is finitely supported, it follows that only finitely many factors 
in (9) are different from 1. 

^In what follows, we also feel free to write the condition determining the set A 
instead of the argument A itself in expressions like Pm{A). 

^If the Qj are integers, then is exactly the number of ways in which 

Nj indistinguishable particles can be placed on an energy level of multiplicity qj, 
and w{{Nj}) is the number of distinct system states corresponding to the sequence 
{Nj} of occupation numbers. 

^The asymptotics (12) with Q = T{d)~^ holds, for example, for the multi- 
plicities Qj = (^'+^"^) of energy levels of the d-dimensional quantum-mechanical 
harmonic oscillator. 

^Most of the results stated below, except for some explicit formulas like (19), 
remain valid if we drop (12) and only require that the estimates (13) be true. 
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has a unique solution b > (which tends to zero as M ^ oo),^ and we 
set 



(15) 



— 1 

In the theorem below, we also need the sum 



(16) 



which tends to infinity together with M. Take a positive function x(x), 
X > 0, tending (arbitrarily slowly) to infinity as x — > oo and set^ 



(17) 



AN\nNfl\{N) ifd>2, 
N^''^hiNx(N) iil<d<2. 



Now we are in a position to state our assertions. 

Theorem 2. Let the numbers qj satisfy (12). Then there exist 
constants > 0, s = 1, 2, . . . , such that the estimates 



(oo 

.7=1 



> A 1 < CsN~ 



1,2,..., 



where Nj, N, and A are defined in (15), (16), and (17), hold for 
an arbitrary M > and an arbitrary sequence {fj}, j = 1,2, ... , of 
complex numbers satisfying the condition sup^- \fj\ < 1 . 

By taking fj = for j < / and fj = 1 for j > / in (18), we obtain 
the following assertion. 

Corollary 1. Under the assumptions of the theorem, for arbitrary 
integer I > 1 one has 



(oo oo \ 

3=1 3=1 ^ 



N~ 



s = l,2,.... 



Essentially, Theorem 2 and Corollary 1 say that for, large M, a ran- 
dom element {Nj} in Xm is well approximated by the Bose-Einstein 
distribution (15) with parameter b. The following closed-form asymp- 
totic expressions relating M, 6, N , and the cumulative distribution 



Indeed, it follows from (13) that the series on the left-hand side in (14) con- 
verges for all 6 > 0; moreover, the sum of this series is easily seen to be a function 
monotone decreasing from oo to as 6 goes from to oo. 

^One should bear in mind that 6, Nj, N, and A are functions of M, even 
though we do not always write out the argument M explicitly. 
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^j>i ^3 as M ^ OO by Proposition 4, (i) in Sec. 4: 
M = h-'^-^QV{d + l)C{d + 1)(1 + o(l)), 
N = r'^Qr(d)C(d)(l + o(l)) and hence 



(19) 



ITT ^ 1 r(r/)(;(f/)(i + o(i)) 



(r(d+i)c(d+i))^ 

oo 

The last formula in (19) reveals the role of the nondimensionalized 
Bose-Einstein distribution function 

rpd—l 

(20) m = -— -. 



It is no surprise that, being appropriately normalized, the random el- 
ements {Nj} e JIm, in a sense, tend as M ^ oo to the function (20). 
More precisely, let us define a sequence of independent random func- 
tions 4>m{x), M = 1, 2, . . . , on the positive real line R+ by setting 

0^(x) = g-V-^AT., :ce[bij-l),bj), j = l,2,..., 

where M and b are related by (14) and {Nj} e Qm is a random element 
of the probabihty space Xm- Then the following assertion holds. 

Corollary 2. The sequence {0m} almost surely *-weakly con- 
verges in the space (C"^(]R+))* of continuous linear Junctionals on the 
space C^(]R+) to the function (20). Namely, for each diff'erentiable 
function f{x), x G bounded together with its first derivative uni- 
formly on R+; one has 

where the angle brackets denote the pairing 



oo 

{u, v) = J 



u{x)v{x) dx. 



The convergence to the limit distribution (20) can also be stated in 
a somewhat different manner, as was done by Vershik [3]. 

Corollary 3 (cf. [3, Theorem 4.4]). For any e > and any closed 
interval [xi,X2\, < Xi < X2 < oo, there exists an Mq such that 



oo 

'm ( sup Q'^b"^ ^ ^j- / ^(^) dx > s) 



< e 



for M > Mq. 
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Remark. Note that the precise information on the convergence 
rate contained in Theorem 2 and Corollary 1 has been lost in Corollar- 
ies 2 and 3. 

2.2. System with a fixed number of particles. Now let us 

consider systems in which two constraints, one on the total energy and 
one on the number of particles, arc given. 

2.2.1. Definition of the system. Let M, > be integers. By Qm,n 
we denote the set of all sequences {Nj} = {Nj}jZ,o of nonnegative 
integers satisfying the conditions^° 



oo oo 



Again, such sequences are finitely supported, and flM,N is finite. We 
take the same multiplicities qj > 0, j = 1, 2, . . . , as in Sec. 2.1 and 
supplement them with some number go ^ 1- Next, we introduce the 
probability space Xm,n = i^M,N, ^m,n,Pm,n), where J='m,n = 2"^'^ 
and the probability Pm,n is defined as follows: 

W{A) 

(21) Pm,n{A) = r, A C ilM,N: 

[}''M,n) 

where, for every A C ^m,n, the weight Ty(^) is given by^^ 

{Nj}eA j=Q ^ 3 / 

2.2.2. Limit distribution and Bose condensation. Let us study the 

behavior of random elements {Nj} G Qm,n as M ^ oo and oo. 
This problem involves two large parameters M and A^ rather than one, 
and it is natural to expect that the answer is more complicated than in 
the case of one large parameter M, considered in Sec. 2.1. It turns out 
that the asymptotic behavior of our sequences {Nj} strongly depend on 
how the rates at which M and A^ tend to infinity are related. Namely, 
let N = N{M) be defined by (14) and (16). RecaU that, by (19), 

N{M) = MH^QHir r(^)y) ^ (1 + o{l)). 

(r(ci+l)C(ci+l))^ 



-'^'^It is merely a matter of convenience that we have decided to start the indexing 

from j — 0, that is, have chosen zero for the ground energy level. Should we wish to 
start from j = 1, it suffices to change the notation as follows: N = N,M = M + N, 
Nj = Nj-i, j = 1,2, . . . . In terms of the variables with tildes, the sums start from 

i = i- 

^^We denote the weights in this section by the capital letter W so as to avoid 
confusion with the weight w introduced in Sec. 2.1 for sequences starting from 
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There are two possible types of asymptotic behavior of random ele- 
ments {Nj} e flM,N as M, N —>■ oo depending on whether N is smaller 
or greater than A'^. 

(i) li N < N{M), then the limit distribution is the Bose-Einstein 
distribution 

(22) ^^ = ^pj^v ^'-O'l'^'---' 

where the parameters (3, /j, > are determined from the system 
of equations 



(23) T-T^ = N, Y-J^^M. 



(ii) li N > N{M), then Bose condensation occurs: the occupation 
numbers Nj with j > 1 in the limit distribution no longer 
depend on and coincide with the numbers (15), while the 
excessive particles, however many, occupy the zero level. 



oo 



(24) Nj^-^, j>l, No^N-J2Nj^N-N. 

Here we do not consider case (i) in detail and refer the reader to 
Theorem 5 in [9], where the probabilities of deviations from the limit 
distribution (22), (23) arc estimated. 

In case (ii), the following theorem holds. 

Theorem 3. Let the numbers qj satisfy (12), and let go ^ 1- Then 
there exist constants Cg > 0, s = 1,2, . . . , such that the estimates 



(25) 



V j=Q / 



s = l,2,.... 



where the Nj are defined in (24) and A is the same as in Theo- 
rem 2, hold for arbitrary M > and N > N{M) and an arbitrary 
sequence {f j} , j = 0,1,2 of complex numbers satisfying the condi- 
tion supji/jl < 1. 

In the same way as in Sec. 2.1.2, we obtain the following corollary. 

Corollary 4. Under the assumptions of the theorem, for arbitrary 
integer I > one has 



M 



3=1 3=1 



> A\<C..N , s = l,2 



In particular, this proves Theorem 1. Indeed, to derive the estimate 
in that theorem from Corollary 4, it suffices to set Z = 1 and notice 



10 



V.P. MASLOV AND V.E. NAZAIKINSKII 



that 



\No-{N -N)\ 



i=i i=i 



Remark. One could also state the counterparts of Corollaries 2 
and 3; we do not dwell upon this. 

3. Proofs 

In this section, the letter C is used to denote various positive con- 
stants independent of A^, b, M, etc. These constants are not assumed 
to be the same in all formulas! If we need to keep track of several 
constants simultaneously, we equip C with subscripts. We also widely 
use the following standard notation: we write / x 51 if the ratio f /g is 
bounded above and below by positive constants; in other words, / and 
g never have opposite signs, f = 0{g), and g = 0{f). 

3.1. Proof of Theorem 2. Instead of (18), it suffices to prove 
that (the modulus sign is removed) 

(26) Pm [Y. - ^^■) >A< CsN~\ s = 1, 2, . . . , 

\j=i / 

where the fj are assumed to be real. Then we obtain (18) for real 
fj (with the constants Cg multiplied by 2) by combining (26) with the 
similar inequality where each /, has the same modulus and the opposite 
sign, and finally reach the case of complex fj using Pythagoras' theorem 
(with further increase in the constants) . So we concentrate on the proof 
of (26). 

By (11), the probability on the left-hand side in (26) is given by 



'M 



w{nM{A)) 



where Qm(A) C is the set of all sequences {Nj} for which 

00 

(27) E/,(iV,-iV,.)> A- 

To prove (26), we will obtain a lower bound for w{D,m) and an upper 
bound for ^((^^(A)). 

3.1.1. A lower bound for w{Qm)- It is not so easy to estimate 
w{flM) directly. Instead, we will estimate the weight w{Q%), where 

C Qm is the subset formed by the sequences for which j^j = 
M (i.e., equality takes place in (8)). This weight obviously does not 
exceed w{Q.m)- First, let us write out an exact formula for w{^\j). Let 



(28) F{z) = J2 M^m) 



M=0 
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be the generating function of the numbers w(O^). 

Lemma 1. The series (28) converges in the disk {\z\ < 1}, and the 
sum is given by the formula 



oo ^ 



zi)<ij ' 



This assertion is well known for integer qj (e.g., see [3] and [10, 
Chap. 1]), and the proof for noninteger qj is also easy.^^ 
Now we can express ^(Q^) by the Cauchy formula 



F{z) dz _ 1 



■ oo ^ 



L (1 _ 



dz 



\z\=r 

where < r < 1. The change of variables z — e~^ yields 
1 



Mj 



2m 



Mi) 



70 



Here the main branch of the logarithm is taken, and the integration 
contour 70 is the circle {Re^ = — Inr} on the cylinder C/27riZ with 
coordinate ^ mod 27ri. To estimate the integral, we use a saddle-point 
argument (cf. [13, Chap. 4]). Let us deform 70 into a contour 7* on 
which 

min max Re $(0 

7 ?e7 

is attained, where the minimum is taken over all contours 7 lying in 
the right half-cylinder {Re^ > 0} and homotopic to 70. We have 

00 ^ 

Re $ (0 = M Re e + 5^ In _. . 

j=i I ^ I 

To find 7*, the following lemma will be of help. 



12t 



Here is the proof. From (28), using (9), (10), and the fact that the sets fl^ 
are disjoint, we obtain 



En 



Nj + - 1 



where the sum is taken over all finitely supported sequences {Nj} of nonnegative 
integers. Since the qj grow polynomially by condition (13), it follows by routine 
estimates that the sum and the product can be interchanged, and we obtain 



oo p oo /7>r , 1 
j=l '-Nj=0 ^ ^ 

by the binomial series formula. 



oo 



LI (1 _ 
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Lemma 2. (i) For fixed Re^ > 0, the maximum o/Re$(^) is at- 
tained on the real axis. 

(ii) The minimum o/Re$(^) on the positive real axis is attained at 
the point ^ — h, where b is the root of Eq. (14). 

Proof, (i) In fact, a stronger assertion is true:^^ if Re^ > 0, then 

^ oo 

(29) Re*(Ree) - Re^O > ^ E ^i^"'""'^!! ' cos(j 1^0)- 

i=i 

All terms on the right-hand side in (29) are nonnegative, and the 
first term (j — 1) is strictly positive unless Im^ = mod 27r. This 
proves (i). 

(ii) For real ^, we have Re$(^) = $(^). Next, 

OO . OO -2 j'^ 

(30) 4>'K) = M-5:-f«L, *"({) = 

j=l j=l V / 

The function $'(0 has the unique positive zero ^ — b (cf. (14)), and this 
zero is the unique point of strict minimum of Re$(^) on the positive 
real axis, because the second derivative is strictly positive for real ^. 
This gives (ii) and completes the proof of Lemma 2. □ 

It follows from Lemma 2 that for 7* we can take the circle {Re^ = 
b} on the cylinder C/27riZ.^^ Hence we write 

IT 

(31) = ±-J e*(« dC^^I e^(-) d^, 



-"^^Indeed, the left-hand side of (29) can be expressed as 
Re $(Re ^ - Re $(0 = ^ Qj In 

Let a; = j Re^ > and y = j Im^. Since Int; > 1 — for w > 1, we have 

|l_e-^-»!/| 1-e-^ VT+S-l . 2e-^ ,^ , 

in- ^ > 1 — 1 ^ = , — , o = - -r l — cosw). 

1-e-^ - ll-e-^-'wl (l-e-=^)2^ 

Next, 



> - — , w = 2e ^(1 - cosy) < (5. 

It remains to note that ui E [0,4], so that 

vTTw-1 ^ yr+w-i uj ui 1 , 

, — > ;= = — r = > tt: = (1 — cosw) 

VT+TJ ~ V5 2VTT^V5 - 10 5 ^ ''^ 

(where 6 G [0, 1]), and we arrive at (29). 

^^Indccd. let 7 bo any contour homotopic to 70. Then 7 necessarily contains a 
point 77 of the real axis, and 

(i) 

maxRe$(0 = $(6) < $(??) < maxRe$(0, 
which shows that 7* is a saddle-point contour. 
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where 

Let us estimate the integral in (31). By construction, we have 
5(0) = 0, and 



oo 

(32) KeS{^) < --J2qje-'^{1 - cos(j(^)) 

by (29). Next, 



oo 



,2 



(33) S'{0) = i^'ib) = 0, S"{0) = -$"(6) = - J] _ 

j=i ^ ^ 

(cf. (30)), and finally. 



IV'J - « 2^ (gbj+uv - 1)3 



Lemma 3. There exist positive constants ci,C2,C3 5«c/i that 

(35) - cife-"^-' < S"{0) < -C2r'^-^ sup|5'"((^)| < csb-"^-^ 

f 

Proof. We use inequalities (13). Then it follows from (33) and 
(34) that 

^ jd+l hj ^ Ad+l bj 

and (since |e''-'+^^'^ -\\>e^^- 1) 

^ j</+2/ ,2/,/ , fejX 

i=i ^ ' 

Now it remains to use Proposition 4, (ii) in Sec. 4. This proves the 
lemma. □ 

We split the integration interval [— tt.tt] in (31) into three zones, 

Di = {|(^| < ^2 = {bvh^^"!^ < \ip\ < d2b}, 

Ds = {52b < \(p\ < tt}. 

Here Si and ^2 are sufficiently small positive constants, independent of 
b, to be chosen later. Let us represent S{(p) by Taylor's formula with 
remainder of order 3, 

S{ip) = ls"{0)ip^ + Rsi^p), where |i?3(<^)| < 
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by the second inequality in (35). In particular, 

(36) |i?3(^)|<^ inD,. 

Take ^3 so small that the right-hand side in (36) is smaller than 7r/4. 
Then 

ce-''^'-'-'^' < Re e^(^) < cge"'^^^"'"'^' in D^ 
for some positive constants C4, . . . , C7. We have 

5i6-i/6 

J e-'^-'-"^" dip = h'l^^^ I e-^y" dy x b'/'+\ 

Di -5i6-i/6 

and hence 

Now consider the zone Z^2- There we have 

ReS{ip) < -(026-'='-^ -C3r'^~-%|)(^2 

< -(C2 - S2C3)b-''-V < -(C2 - S2C3)b-''/'Sl 

and if we take S2 small enough that cs — C2 — S2CS > 0, then the integral 
over D2 decays exponentially (at the rate of e"'^*'' ''^^) as 6 — 0. 

Finally, consider the zone D^. Here we use inequality (32). Since 
all terms in the series on the right-hand side in (32) are nonnegative, 
we can drop some terms and write 

ReS{ip)<~ Yl qje-'^0--cos{jip)) 

xi<bj<X2 

with some positive Xi and X2 to be chosen later. By (13), qje^''^ > 
096"'^+^ for xi < bj < X2, where the constant Cg depends on xi and X2, 
and so 

ReS{ip)<-^^ E (l-cos(j^)). 

Xl<bj<X2 

Using [12, 1.341.3], we obtain 
Yl " cos(j99)) 

Xl<bj<X2 

(37) _ . COS ((ji +i2)¥'/2) sin ((j2-ii + 1)^2) 

— J2 — Jl + J- H r-7 TTTT 

sin((/?/2) 

> b-\x2 - xi) - 1 - \sm{<p/2)\-\ 

where ji and j2 are the first and the last integer, respectively, in the 
interval (xi/6, X2/6). Since cp e D3, we have |sin((/7/2)| > sm{62b/2) > 
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82b /ti, and so the sum (37) is greater than Cio&~^ provided that we take 
a sufficiently large X2 — Xi. Thus, we see that 

ReS{ip) <ciih''^ inDa, 

so that the integral over also decays exponentially as 6 — > 0. Now 
we summarize the preceding and see that we have proved the following 
assertion: 

Proposition 1. There exists a constant C > such that 
where 

oo ^ 

(38) S{M) ^bM + J2 Qj In YZr^ 

and b is related to M by formula (14). 

This is the desired lower bound for w{Qm)- 

3.1.2. An upper bound for w{Qm{^))- It follows from (8) and (27) 
that if {Nj} e fiM(A) and c > 0, then 



(39) 



y OO \ OO 

b(M-J2jNA+cJ2 fj (Nj - F,) - cA > 0. 

V 7 = 1 ^ 7 = 1 



Accordingly, the exponential of the left-hand side of (39) is greater 
than 1 for any {Nj} G Qm{^)', of course, it is positive for any (finitely 
supported) {A'^}. Hence 



(40) 

{Nj} 



— e 



En 



Nj + Qj - 1 \ (-bj+cfj)Nj 

AT. 



where the sum extends over all finitely supported sequences {A^} of 
nonnegative integers. Note that if c < 6/2, then 



e 



-bj+cfj ^ (^(,-b/'^y 



(recall that \fj\ < 1), and hence the series on the right-hand side in (40) 
is dominated by the series (28) with z — e~^^'^ < 1, which is convergent 
by Lemma 1. Arguing as in the proof of Lemma 1, we interchange the 
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sum and the product and use the binomial formula to obtain 



w ^ — . — 

'-Ni=0 

oo 



\-M —n V 3 / 



oo 



L (1 — e-^3+cijYi 

expifeM - cA + y^g, In r^^^ 1 

By Taylor's formula with remainder, 



where 9j G [0, 1]. Since \fj\ < 1 and c e [0, 6/2], we obtain, transposing 
the term cfj{e''^ — 

cfj 1 , c2 e'^iJ+y^^ 



In r-— ? 1—— < In + 



1 _ (,-bj+cfj gfej _ 1 - 1 _ g-6i 2 (e^(J-V2) _ 1)5 

1 c2 e*e*^/2 

< In + 



1-e-^J' 2 (e''^/2-l)2' 
and thus we arrive at the following assertion. 
Proposition 2. One has the estimate 

w{nM{/^)) < e'5We-^^+'='^, c e [0, 6/2], 
where S{M) is given by (38) and 

2 (e6j/2_i)2- 

3.1.3. Completion of the proof . By combining Propositions 1 and 2, 
we obtain 

Pm fj{Nj - Nj) >a\< Ch-^-'l-'e-'^^^'^''', c e [0, h/2]. 

Now we will prove that there exists a c e [0, 6/2] such that 
(41) e-cA+c^^ < s = 1, 2, . . . . 

Then the assertion of the theorem readily follows, since 6 x 
We consider three cases. 

1. d > 2. In this case, K < Ch~''' by Proposition 4, (ii) in Sec. 4, 
and soK < CN. Set 
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Since N x b~^, it follows that c < b/2 for sufficiently large M. Next, 
and (41) holds. 

2. d — 2. In this case, K < Cb^'^\lnb\ (Proposition 4, (ii)) and 
hence K < CNlaN. Set c = b/2. Then 

cA - c'^K > ^VTlnNx(N) - CjNlnN. 

The first term is of the order of lnNx{N), and the subtrahend is of 
the smaller order of In , so that (41) again holds. 

3. 1 < d < 2. In this case, K < Cb~'^ (Proposition 4, (ii)), and 

hence K < Cn'^''^ . We again take c = b/2. Then 

cA - c^K > p^'^'lnNxiN) - C^^N^^" . 

The first term is of the order of In A'^x(A^), and the subtrahend is 0(1), 
so that (41) again holds. 

The proof of Theorem 2 is complete. □ 

3.2. Proof of Corollary 2. The proof of this corollary is based 
on the following lemma. 

Lemma 4. One has 

oo 

(42) (0 - 0M, /) = h-'Q-^ ^(iV,- - Nj)f{bj) + R,{M) + R^{M), 
where Ri{M) is a random variable on such that 

(43) \Ri (M) I < Cb'^+^ - £ Nj) 

and R2{M) is a {nonrandom) function such that 

R2{M) ^0 as M ^ oo. 

Proof. Since the function f{x) is bounded together with the first 
derivative and 0(x) is given by (20) with d > 1, it follows from the 
Euler-Maclaurin formula (Proposition 3) that 

/oo 
0(x)/(x)da; = bJ2<f^{bj)f{bj)+r{M), 

where r(M) ^ as M ^ oo. Next, by (12), 
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where ^ as j — > oo. We claim that 

oo 

b'^Q-^ J2 Njejfibj) ^0 as M ^ oo. 
Indeed, since / is uniformly bounded, we have 

oo / oo N 

b'Q-'Y.^jejfibj) <cb''lj2Nj\ej\+ 

j=l ^j=l j=k+l ^ 

where k is arbitrary. Take k so large that \Qj\ < e for j > k. Then the 
second sum does not exceed eN. For fixed k, the first sum does not 
exceed Cib~^ (where the constant depends on k) in view of formula (15) 
for Nj. Thus, the right-hand side does not exceed CCib'^'^ + Cb'^'Ne. 
Recall that N x b'"^] hence the second term can be made as small as 
desired by an appropriate choice of e, and then we take M large enough 
(i.e., b small enough) to ensure that the first term is also small. We 
conclude that 

oo 

(44) (0, /) = b'Q-'Y.N^fibj) + ri(M), 

3=1 

where ri(M) — * as M ^ oo. 

Now let us study {(f)M,f)- We represent / in the form f[x) = 
fi{^) + f2{x), where fi{x) is the step function taking the value f{bj) 
on the interval [b{j — l),bj) and f2{x) — f{x) — fi{x) satisfies the 
estimate 1/2(1^) | < Cb, since the derivative f'{x) is uniformly bounded. 
Then, in view of the definition of 4>m{x), we have 

00 

(0M, /) = (0M, /l) + (0M, /2) = b^Q-' Yl ^^•^(^^■) + ^(^)' 

where 



00 

/oo 00 
(I)m(x) dx = Cb'^+^Q-^ Nj = Cb'^+^Q-^J2^. 
3=1 

+ Cb'^+^Q-^ (n - = 
V ,=1 / 



3=1 3=1 
r2(M) +r3(M). 



Now we can set -Ri(M) = r^iM) and R2{M) = ri(M) + r2(M). The 
proof of the lemma is complete. □ 

We assume without loss of generality that |/(a;)| is bounded by 1. 
It follows from Theorem 2 and Lemma 4 that for each £ > there 
exists an Mq = Mo{e) such that 

(45) Pm(|(0-0m,/)|>£)<CM-2 forM>Mo. 
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Indeed, the third term on the right-hand side in (42) is necessarily less 
than e/3 in modulus for sufficiently large M. Let us study the ffist 
term. It is necessarily less than s/3 in modulus provided that 

oo 

5]/(6j)(7V,-]v,) <C,N, 

where Cq is some constant (depending on Q and on the constants in 
the relation x 6^'^. If M is large enough that CqN > A (which can 
always be achieved, because A grows slower than N as M ^ oo), we 
can apply Theorem 2 with appropriate s and conclude that 



M 



> £/3 < CM-\ 



i=i 

The same argument applies to the second term on the right-hand side 
in (42), and we arrive at (45) (with some new C and Mo(£)). 

Consider the product nM=i '^m of the probability spaces Xm- The 
probability measure on this product will be denoted by P. Consider 
the event 

A^e = {\{<p-(t>Mj)\<e for all M>m}. 
It follows from (45) that 

oo 
M=m 

for m > Mo{e). The series on the right-hand side converges, and so 

P{Ame) — 1 as m ^ OO. 

Next, for the probability of the desired convergence (0m, /) — ^ (0> /) 
we have the expression 



P((0 



J) ^ 0) = P(V£ > 03m: A^,) = p(n U • 



Note that the sets Ay^y^ are nested, 

Arae C A^a'e for m < m\ 

and that the sets 

^£ = Ams 

are nested as well, 

B, C for e < e'. 

It follows that 

P{B,) = hm P{A^,) = 1, P((0 - (Pm, /) ^ 0) = lim P{B,) = 1. 

m— »oo £— ►0 

The proof of Corollary 2 is complete. □ 
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3.3. Proof of Corollary 3. Using Proposition 4 and arguing as 
in the proof of (44), we obtain 

oo 

J <j){x) dx = b'^Q-^ + ^(^)' 

J. j>x/b 

where r(M) — > as M — > oo. Hence it suffices to prove that 




where C is a constant depending on Q and the constants in the relation 

X b~'^. For sufficiently large M, we have eCN > A, and Theorem 2 
can be used to estimate the probability for each fixed x. This probabil- 
ity is less that CgN * for all s = 0, 1, 2, ... . The supremum is actually 
taken over the finite set of integer points on the interval [xi/b,X2/b], 
and the number of these points is of the order of 1/6, i.e., of the order 

of A^ ' . Thus, passing from individual points to the supremum over 

0{N^^'^) points, we make the estimate of the probability slightly worse 
(by the factor equal to the number of these points); i.e., the probability 

estimate becomes CgN ' with some new constants Cs- However, 
this does not matter, and we take s = 1, which provides the desired 
estimate: 




< e 



for sufficiently large M (and hence A^). 
The proof of Corollary 3 is complete. 



3.4. Proof of Theorem 3. Just as in the proof of Theorem 2, it 
suffices to show that 

(46) Pm,n (f^ fAN, - iV,) >a\< C,N-\ s = 1,2,..., 

i.e, drop the modulus sign. Next, we will assume without loss of gen- 
erahty that /o = 0. Indeed, 

oo oo 

i=o j=o 

and hence 

oo oo 

j=0 j=l 
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The new numbers /j = fj — fi arc bounded in absolute value by 2 rather 
than 1, but we can pass to the numbers fj = /j/2 and simultaneously 
divide the function x{^) occurring in definition of A by 2. 

By definition (21), the probability on the left-hand side in (46) has 
the form 



(47) PM,iv(f;/,(iV,-7V.)>A^ 



W{Qm,n{A)) 

w{nM,N) ' 

where Qm,n{^) C Om.at is the set of all sequences {Nj}'jLQ e flM,N 
such that 

oo 

(48) ^/,(iV,-iV,)> A. 

Consider the mapping 

r : ^M,N — ^ 

that takes each sequence {Nj} = {Nj}'^^ £ ^m,n to the sequence 
{NjY = {Nj}'jLi G Qm obtained by throwing away the first ele- 
ment Nq. This mapping is one-to-one onto its range, because Nq = 
N — Nj is uniquely determined by {Nj}'. Since /o = 0, it follows 
that 

r((]M,7v(A)) C 1^m(A). 

(Inequality (48) does not involve Nq.) 
For each {Nj} e ^m,N: we have 

w(W})=(''°+^;"')»(W}'). 

Since, by assumption, gfo > 1, we see that the binomial IS a 

monotone increasing function of A^o; by Stirling's formula. 



Now we can write 



4^C°T')] .5: 



{Nj}enM,N{A) 



< CN'i°-^w{nM{A)). 
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On the other hand, 



{Nj}enM,N ^ / 



{NjjeQM, 

We have N > N. Let us consider two cases, > + A and 

N < N <N + 

1. Let > + A. The set \ t{Q,m,n) consists of all sequences 
{Nj} e flM for which 

oo oo 

j=l j=l 

Hence 

— — '■ = Pm{w{Qm) - w(t(Qm,7v))) < CsN , 

w{Um) 

s = 1, 2, .... by Corollary 1 (with / = 1). It follows that, for sufficiently 
large M and A^ > N{M) + A, 



2 I^HM.iv y A'q 
and 

(49) PM,iv(E/,(A^, -A^,)> a) 
\j=i / 

If A^ < M, then the minimum is attained at A^o = and is equal to 1. 
If A^ > M, then the minimum is attained ed, Nq — N — M, and we 
have, by Stirling's formula, 



In any case, we obtain 



Pm,n (f2 MN, - N,) >a\< CsN 
V,=i / 
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(with some new constants Cg). This proves the theorem for this case. 

2. Let < < iV + A. Here we need a finer argument. 

We carry out the same computations for the upper bound of the 
numerator, but to estimate the denominator, we write 

p^(^)M,.)>('^-J+|-') E n^m}) 

fN-N + qo-l\ ^ 
where fl'^ C flu is the subset of sequences satisfying the conditions 

oo oo 

J] AT,- = TV, j;jW, = M. 

j=i i=i 

It follows from Lemma 6 and formula (40) in [9] that, for some mo, 

where S{M) is defined in (38). By combining this with the estimate 
for w{riM{^)) obtained in the proof of Theorem 2 (Proposition 2), we 
arrive at the desired estimate. 

The proof of Theorem 3 is complete. 



4. Appendix 

We shall use the following two versions of the Euler-Maclaurin for- 
mula. 

Proposition 3. Let f{x) be a continuously dijferentiable function 
on (0, oo) absolutely integrable together with f'{x) at infinity. Then the 
series Yl'jLi fU) converges absolutely, and its sum can be computed by 
the formula 



^ CO 

oo „ 

5^/(j)= / f{x)dx + R,, 



(50) 

where the remainder Ri satisfies the estimate 

oo 

\Ri\ < J \f{x)\dx. 

1 

If in addition, f'{x) is absolutely integrable at zero, then 

oo p 

(51) E^(^')= / f{^)dx + R2, 

7 = 1 { 
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where the remainder R2 satisfies the estimate 



00 

m < J \f'{x)\dx. 



Proof. The obtain the desired error estimates, one can use the 
Stieltjes integral representations 

00 00 

•5=^ 1 

where [x] is the integer part of x, and then integrate by parts to estimate 
the remainders. □ 

Let us use Proposition 3 to carry out some computations needed in 
the main text. 

Proposition 4. The following asymptotic formulas hold as b ^ 

+0. 

(i) Let s > 0. Then 
00 .g 

E = + i)c(. + i)(i + 0(1)) 

and, more generally, 



00 



^ e^-? - 1 ev -1 ^ ' 



(ii) Let d > 1. Then 



Z-^(ef^_l)2 J 

00 



Cb-'^ ifd>2, 
T.7^EJJ^2<\Cb-'\lnb\ zfd^2, 
J=i ^ ' [C6-2 ifl<d<2 

Proof, (i) If s > 1, then, by (51), 

00 

f dx L-s-i /" 



00 00 



J=i 

= r*-^r(s + i)C(s + i) + i?2 



ev-l ' 
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(we have used [12, 3.411.1]), where 







and we have the desired estimate, because the integral converges. If 
< s < 1, then this computation does not work, but we can use (50) 
and write 



oo oo 



1 h 

oo b 



= 6 



—s- 



J ev-l J ey-1 







where 



b 

and so < Cb~^ for s G (0, 1) and < Cb-^\lnb\ for s = 1. Now 
note that 

' ^ < 

ey -1 ~ 



and we obtain the desired estimate. 

The estimate for the sum starting form j — I can be obtained in a 
similar way; one only need to shift the summation index by Z — 1. 

(ii) The proof goes by the same scheme as for (i). We omit the 
details. □ 

Remark. The integrals in (ii) can also be expressed via special 
functions, but we do not need these expressions. 

References 

[1] L. D. Landau and I. M. Lifshits, Theoretical Physics, Vol. 4: Statistical Physics 

(GITTL, Moscow-Leningrad, 1951) [in Russian]. 
[2] I. A. Kvasnikov, Statistical Physics. Vol. 2: Theory of Equilibrium Systems 

(URSS, Moscow, 2002) [in Russian]. 
[3] A. M. Vershik, "Statistical mechanics of combinatorial partitions, and their 

hmit shapes," Funktsional. Anal, i Prilozhen. 30 (2), 19-39 (1996) [Functional 

Anal. Appl. 30 (2), 95-105 (1996)]. 
[4] A. M. Vershik, "Limit distribution of the energy of a quantum ideal gas from 

the viewpoint of the theory of partitions of natural numbers," Uspekhi Mat. 

Nauk 52 (2), 139-146 (1997) [Russian Math. Surveys 52 (2), 379-386 (1997)]. 
[5] V. P. Maslov and V. E. Nazaikinskii, "On the distribution of integer random 

variables related by a certain linear inequality, I," Mat. Zametki 83 (2), 232- 

263 (2008) [Math. Notes 83 (2), 211-237 (2008)]. 



26 



V.P. MASLOV AND V.E. NAZAIKINSKII 



[6] V. P. Maslov and V. E. Nazaikinskii, "On the distribution of integer random 

variables related by a certain linear inequality, II," Mat. Zametki 83 (3), 381- 

401 (2008) [Math. Notes 83 (3), 345-363 (2008)]. 
[7] V. P. Maslov and V. E. Nazaikinskii, "On the distribution of integer random 

variables related by a certain linear inequality. III," Mat. Zametki 83 (6), 

880-898 (2008) [Math. Notes 83 (6), 804-820 (2008)]. 
[8] V. P. Maslov and V. E. Nazaikinskii, "On the distribution of integer random 

variables related by two linear inequalities, I," Mat. Zametki 83 (4), 559-580 

(2008) [Math. Notes 83 (4), 512-529 (2008)]. 
[9] V. P. Maslov and V. E. Nazaikinskii, "On the distribution of integer random 

variables satisfying two linear relations," Mat. Zametki 84 (1), 69-98 (2008) 

[Math. Notes 84 (1), 73-99 (2008)]. 
[10] G. E. Andrews, The Theory of Partitions (Addison- Wesley, Reading, 1976). 
[11] A. A. Karatsuba Fundamentals of Analytic Number Theory (URSS, Moscow, 

2004). 

[12] I. S. Gradshtein and I. M. Ryzhik, Table of Integrals, Sums, Series, and Prod- 
ucts (GIFML, Moscow, 1963) [in Russian]. 

[13] M. V. Fedoryuk, The Saddle-Point Method (Nauka, Moscow, 1977) [in Rus- 
sian]. 

Faculty of Physics, Moscow State University 

Institute for Problems in Mechanics, Russian Academy of Sciences 

E-mail address: v.p.inaslov@inail.ru 

Institute for Problems in Mechanics, Russian Academy of Sciences 
E-mail address: nazaikinskii@yandex.ru 



